Nutrients 2014, 6, 2987-2999; doi:10.3390/nu6082987 



OPEN ACCESS 



nutrients 

ISSN 2072-6643 

WW w . mdp i . c om/j oumal/ nutrients 

Article 

Association between Intake of Sugar-Sweetened Beverages and 
Circulating 25-Hydroxyvitamin D Concentration among 
Premenopausal Women 

Caroline S. Duchaine ^'^ and Caroline Diorio 

^ Cancer Research Center, Department of Social and Preventive Medicine, Laval University, Quebec 
City, Quebec GIV 0A6, Canada; E-Mail: caroline.duchaine.l@ulaval.ca 
Oncology Unit, CHU de Quebec Research Center, Saint- Sacrement Hospital, Quebec, 
QCG1S4L8, Canada 

Deschenes-Fabia Center for Breast Diseases, Saint-Sacrement Hospital, Quebec, QC GIS 4L8, Canada 

* Author to whom correspondence should be addressed; E-Mail: Caroline. diorio@uresp.ulaval.ca; 
Tel.: +1-418-682-7511 (ext. 4726); Fax: +1-418-682-7949. 

Received: 5 May 2014; in revised form: 11 July 2014 / Accepted: 21July 2014 / 
Published: 28 July 2014 

Abstract: Intake of sugar-sweetened beverages has increased in North America and seems 
to have several adverse health effects possibly through decreased circulating 
25-hydroxyvitamin D (25(OH)D) concentrations. The aim of this cross-sectional study was 
to evaluate the association between sugar-sweetened beverages intake and 25(OH)D 
concentrations among premenopausal women. Intake of sugar-sweetened beverages 
including colas, other carbonated beverages and sweet fruit drinks was assessed using a 
validated food frequency questionnaire among 741 premenopausal women. Plasma 
concentrations of 25(OH)D were quantified by radioimmunoassay. The association 
between sugar-sweetened beverages intake and 25(OH)D concentrations was evaluated 
using multivariate generalized linear models and Spearman correlations. A higher intake of 
colas was associated with lower mean 25(OH)D levels (67.0, 63.7, 64.7 and 58.5 nmol/L 
for never, <1, 1-3 and >3 servings/week, respectively; r = -0.11 {p = 0.004)). A 
correlation was observed between intake of other carbonated beverages and 25(OH)D 
concentrations but was not statistically significant (r = -0.06 (p = 0.10)). No association 
was observed between intake of sweet fruit drinks and 25(OH)D concentrations. This study 
suggests that high intake of colas may decrease 25(OH)D levels in premenopausal women. 
Considering the high consumption of these drinks in the general population and the possible 
consequences of vitamin D deficiency on health, this finding needs further investigation. 
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1. Introduction 

Consumption of sugar-sweetened beverages has increased considerably in the past 30 years and has 
become one of the main types of beverages in North America [1]. In the United States, a citizen 
consumes on average 12 servings of soda per week of which 58% contain sugar [2]. In Canada, soft 
drinks have been the leading beverage choice between 1980 and 2000 [1] and they remain so today [3]. 
This excessive intake of sugar-sweetened beverages may have several adverse effects on human health 
such as diabetes, hypertension [2], obesity, cardiovascular diseases [4], low bone mineral density [5], 
hypocalcaemia [6,7] and some types of cancer [8-11]. 

Interestingly, the link between insufficient circulating 25-hydroxyvitamin D (25(OH)D) 
concentrations and similar health outcomes was established with low bone mineral density [12], and 
proposed with diabetes [13], metabolic syndrome [14], hypertension [15] and cancer [16-21]. The 
25(OH)D is the principal circulating vitamin D metabolite and a well-established biomarker for total 
vitamin D exposure from food, supplements and endogenous synthesis [22]. The Institute of Medicine 
report on dietary reference intake for calcium and vitamin D recommends 25(OH)D concentrations 
of 50 nmol/L or more for good health and establishes the deficiency limit at 30 nmol/L [12]. In 
Canada, 25.7% of the population had plasma concentrations of 25(OH)D below 50 nmol/L in 201 1 [23]. 
Thus, an inter-individual fiuctuation in 25(OH)D concentrations may modulate the risk of some 
diseases, and this fluctuation could be due to genetic or environmental/lifestyle factors including 
nutrition [22]. 

In fact, some studies found that intake of sugar-sweetened beverages could have an impact on 
circulating concentrations of 25(OH)D. Sugar-sweetened beverages, particularly colas, have been 
shown to decrease concentrations of plasma 25(OH)D in rats [24]. Following a two-month diet, rats 
consuming cola beverages had significantly lower plasma 25(OH)D concentrations than rats 
consuming water {p < 0.001). Furthermore, one cross-sectional study among children has shown an 
inverse association of intake of soda {p < 0.001) and juice {p = 0.009) with serum concentration 
of 25(OH)D [25]. Some of the nutrients contained in sugar-sweetened beverages, like fi^lctose, 
caffeine and phosphoric acid, have already been proposed to explain the link of such beverages with 
circulating concentrations of 25(OH)D [24,26-28]. To further explore this link, we have evaluated the 
association of intake of sugar-sweetened beverages (colas, other carbonated beverages and sweet fruit 
drinks separately or combined) with plasma concentrations of 25(OH)D in a population of 
741 premenopausal women. 
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2. Experimental Section 

2.1. Study Population and Recruitment Procedure 

Study population and recruitment procedure have been described elsewhere [29,30]. Briefly, 
premenopausal women were recruited between February and December 2001 at a private radiology 
clinic among those who received a screening mammogram. Exclusion criteria were personal history of 
cancer or breast surgery, having used hormonal derivatives within three months prior to 
mammography, having used tamoxifen or raloxifene, being pregnant, having diabetes mellitus, 
dwarfism/acromegaly, thyroid, adrenal, or hepatic disease. Among the 777 premenopausal women 
recruited, 741 provided a written informed consent to use their blood sample for assays other than 
those planned at recruitment. This study was approved by the Research Ethics Review Board-Hopital 
du Saint-Sacrement du CHU de Quebec. 

2.2. Data Collection 

Blood sample, height, weight, and waist and hip circumferences were collected at recruitment by a 
trained nurse. Information on smoking status, alcohol intake and education were collected during a 
phone interview performed by trained interviewers. Physical activity was assessed using the Nurses' 
Health Study II Activity and Inactivity Questionnaire [31] and then expressed as metabolic equivalent 
hours per week (MET-h/week) [32]. 

2.S. Diet Assessment 

Information on diet was documented with a 161 item validated [33,34], self-administered 
semi-quantitative food frequency questionnaire (97GP copyrighted at Harvard University, Boston, MA, 
USA). Nutrient content of food was assessed at Harvard University, using U.S. Department of 
Agriculture sources, supplemented with data from food manufacturers and personal communications 
with laboratories. A standard serving size was specified for each item and its frequency of 
consumption over the previous year, ranging from "never" to "more than three servings per day", was 
reported. From this food frequency questionnaire, four items were classified as sugar-sweetened 
beverages: colas with sugar containing caffeine, colas with sugar but caffeine-free, carbonated 
beverages with sugar other than colas and sweet fruit drinks. 

2.4. Assessment of Plasma 25(OH)D 

At the time of collection, samples of blood were rapidly aliquoted and stored at -80°C for 
subsequent analyses. Plasma 25(OH)D concentrations (nmol/L) were quantified by radioimmunoassay 
after acetonitrile extraction (DiaSorin Inc, Stillwater, MN, USA) between November 2005 and 
January 2006. The results met the performance target set by the International 25-Hydroxyvitamin D 
External Quality Assessment Scheme Advisory Panel in 2004—2005. Four blinded duplicates on 
average were made for each of the 24 batches and the intra-batch and between-batch coefficients of 
variation were 7.3% and 8.8%, respectively. 
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2.5. Statistical Analyses 

The frequency of consumption of each sugar-sweetened beverage was fransformed in number of 
servings per week using the middle of each category when necessary. For example, when the response 
category was 1-3 servings per week, we assigned two servings per week to this data. The number of 
servings per week of each item was categorized as never, <1, 1-3 and >3 servings/week or to ensure at 
least 5% of the population in each group. Since very few women consumed one serving of 
caffeine-free cola or more per week {n = 11), this item was combined with regular cola intake. The 
sum of servings per week of all sugar-sweetened beverages was also calculated and categorized 
according to the same criteria. Crude and adjusted mean concentrations of 25(OH)D for each category 
of servings per week were estimated using generalized linear models, and the /j-value for each estimate 
compared to the reference category is provided. Spearman coefficients were used to assess the 
correlation between intakes of each sugar-sweetened beverage item and 25(OH)D concentrations as 
continuous variables. Determinants of 25(OH)D concentrations identified a priori [35-37] and 
included in the adjusted models were age (year) at recruitment, body mass index (BMI) (kg/m"), 
waist-to-hip ratio, total vitamin D (Ul/day) and total calcium (mg/day) intakes from food and 
supplements in the past year, total caloric (kcal/day) intake in the past year, season at blood sample 
collection and leisure-time physical activity (MET-h/week) in the past year (as proxy variables of sun 
exposure), alcohol intake (servings/week) in the past year, education (highest completed degree: 
primary, secondary, college, university), and smoking status (non, former or current smoker). In order 
to account for the possible displacement of milk in the diet by sugar-sweetened beverages. Spearman 
correlations fiuther adjusted for milk intake were also estimated. Women with missing values were 
excluded from analyses. All tests were 2-sided and a/» < 0.05 was considered statistically significant. 
All statistical analyses were carried out using SAS 9.3 (SAS Institute, Inc., Gary, NC, USA) 
softw^are system. 

3. Results 

Characteristics of the study population are shown in Table 1 . The mean age at recruitment was 
46.8 years and the women presented a BMI of 25.2 kg/m on average. The mean plasma concentration 
of 25(OH)D was 65.0 nmol/L, 23.2% {n = 111) of the women had 25(OH)D concenfrations below 
the 50 nmol/L recommended and 1.75% {n = 13) were below the deficiency level of 30 nmol/L. 
Regarding calcium intake, the mean concentration in this study population (974.6 mg/day) was below 
the dietary reference intake recommended by the Institute of Medicine report (1000-1200 mg/day for 
women aged between 30 and 70 years). Most of the women (43.2%) in this population were 
non-consumers of colas (regular or caffeine free) or consumed less than one serving per week of this 
type of beverage (36.7%). Similar patterns were observed for other carbonated beverages and sweet 
fruit drinks with 44.2% and 50.7% for non-consumers and 44.3% and 35.6% for less than one serving 
per week, respectively. 
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Table 1, Characteristics of the study population. 



Premenopausal women (« = 741) ' 

Age (year), mean (SD) 46.8(4.6) 

Body mass index (kg/m^), mean (SD) 25.2 (4.5) 

Waist-to-hip ratio, mean (SD) 0.78 (0.06) 

Plasma 25(OH)D concentrations, (nmol/L), mean (SD) 65.0 (19.6) 

Vitamin D intake ^ (lU/day), mean (SD) 284.3 (23 1 .6) 

Calcium intake ^ (mg/day), mean (SD) 974.6 (435.0) 

Milk intake ^ (servings/week), mean (SD) 5.5 (5.9) 

Caloric intake (kcal/day), mean (SD) 1 906 (5 14) 

Physical activity (MET-h/week), mean (SD) 27.2 (22.2) 

Alcohol (servings/week), mean (SD) 3 .4 (3 .8) 
Total colas (servings/week), n (%) 

Never 319(43.2) 

<1 271 (36.7) 

1-3 103 (13.9) 

>3 46 (6.2) 

Other sweet carbonated beverages (servings/week), n (%) 

Never 327 (44.2) 

<1 328 (44.3) 

>1 85(11.5) 

Sweet fruit drinks (servings/week), n (%) 

Never 375 (50.7) 

<1 263 (35.6) 

>1 101 (13.7) 

Total of all sugar-sweetened beverages ^ (servings/week), n (%) 

Never 130 (17.7) 

<1 304(41.3) 

1-3 154 (20.9) 

>3 148 (20.1) 

Season at blood sampling, n (%) 

Winter 90 (12.2) 

Spring 249 (33.6) 

Summer 184(24.8) 

Fall 218 (29.4) 

Education, n (%) 

Less than a high school degree 48 (6.5) 

High school degree 235 (3 1 .7) 

College degree 197 (26.6) 

University degree 261 (35.2) 
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Table 1. Cont. 



Premenopausal women (n = 741) 



Smoking status, n 
Nonsmoker 
Ex-smoker 
Current smoker 



337 (45.5) 
294 (39.7) 
110(14.8) 



Abbreviations: 25(OH)D = 25-hydroxyvitamin D, lU/day = International unit per day, MET-h/week = 
metabolic equivalent hours per week. ' Missing values for caffeine free colas (« = 2), other carbonated 
beverages (« = 1), sweet fruit drinks (« = 2). ^ Total intakes from food and supplements. ^ Including skim 
milk, 1% and 2% milk and full cream milk. '* Including regular colas and caffeine free colas. ^ Including 
regular colas, caffeine free colas, other carbonated beverages and sweet fruit drinks. 

Table 2 shows relations of intake of sugar-sweetened beverages and plasma 25(OH)D 
concentrations. In adjusted models, a higher intake of colas was associated with lower mean 
concentrations of plasma 25(OH)D (67.0, 63.7, 64.7 and 58.5 nmol/L for never, <1, 1-3 and >3 
servings/week, respectively; r = -0.11 {p = 0.004)). Consumption of carbonated beverages other than 
colas was negatively correlated with the concentrations of 25(OH)D (r = -0.06), but the association 
was not statistically significant (/? = 0.10). No association was observed of intake of sweet fruit drinks 
or total intake of sugar-sweetened beverages with plasma concentrations of 25(OH)D (r = -0.03 
(p = 0.48) for fruit drinks and r = -0.06 (p = 0.09) for total of beverages). Further adjustment for milk 
consumption did not materially alter our results. 



Table 2. Relations of intake 
25(OH)D concentrations. 



of sugar-sweetened beverages with plasma 



Sugar-sweetened beverages 
(Servings/week) 


n 


Crude 


Adjusted * 


25(OH)D (nmol/L) 
mean (95% CI), 


25(OH)D (nmol/L) 
mean (95% CI), 


Total of all colas ^ 








Never 


319 


67.5 (65.3-69.6), ref. 


67.0 (65.2-68.9), ref. 


<1 


271 


63.6 (61.3-65.9), 0.02 


63.7(61.8-65.7), 0.02 


1-3 


103 


66.0 (62.3-65.9), 0.51 


64.7 (61.4-67.9), 0.21 


>3 


46 


54.4 (48.9-60.0), <0.0001 


58.5 (53.4-63.6), 0.003 


Spearman r (p) 




-0.13 (0.0004) 


-0.11 (0.004) 


Spearman r(p)'^ 






-0.11 (0.005) 


Other sweet carbonated beverages 








Never 


327 


66.4 (64.3-68.5), ref. 


66.1 (64.3-67.9), ref. 


<1 


328 


64.1 (62.0-66.3), 0.14 


64.1 (62.3-65.8), 0.12 


>1 


85 


62.6(58.4-66.7), 0.11 


63.7 (60.0-67.4), 0.26 


Spearman r (p) 




-0.08 (0.03) 


-0.06 (0.10) 


Spearman r(p)'^ 






-0.06 (0.13) 
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Table 2. Cont. 



Sugar-sweetened beverages 
(Servings/week) 




Crude 


Adjusted ^ 


n 


25(OH)D (nmol/L) 
mean (95% CI),p ^ 


25(OH)D (nmol/L) 
mean (95% CI),p ^ 


Sweet fruit drinks 








Never 


375 


65.9 (63.9-67.8), ref. 


65.4 (63.6-67.0), ref. 


<1 


263 


63.2 (60.8-65.5), 0.08 


64.0 (62.0-66.0), 0.31 


>1 


101 


66.2 (62.4-70.0), 0.87 


65.6 (62.3-68.8), 0.91 


Spearman r (p) 




-0.04 (0.33) 


-0.03 (0.48) 


opCailllall A \P} 






-0.03 (0.47) 


Totol of oil siif^or-swpptPTipd hevproi^ps ^ 








Never 


130 


66.5 (63.1-69.9), ref. 


65.5 (62.6-68.4), ref 


<1 


304 


65.1 (62.9-67.3), 0.49 


65.7 (63.8-67.6), 0.91 


1-3 


154 


64.6(61.5-67.7), 0.42 


63.8(61.1-66.4), 0.38 


>3 


148 


64.0 (60.8-67.1), 0.28 


64.3 (61.4-67.1), 0.56 


Spearman r (p) 




-0.09 (0.02) 


-0.06 (0.09) 


Spearman r(p)'^ 






-0.07 (0.08) 



Abbreviations: 25(OH)D = 25-hydroxyvitamin D, ref. = reference group. ' Adjusted for age, body mass index, 
waist-to-hip ratio, physical activity, alcohol intake, total vitamin D, calcium and caloric intakes, season at 
blood sample, smoking status and education. ^ p values obtained from generalized linear models with / test 
comparison between means of 25(OH)D concentrations for each category of sugar-sweetened beverages 
intake and the reference group. ^ Including regular colas and caffeine free colas. " Including skim milk, 1% 
and 2% milk and frill cream milk. ^ Including regular colas, caffeine free colas, other carbonated beverages 
and sweet fruit drinks. 



Table 3. Partial correlations between determinants of 25(OH)D concentration and plasma 
25(OH)D concentrations. 



Determinant of 25(OH)D concentration 




Spearman correlations 




r' 




r' 




Age (year) 


-0.05 


0.18 


-0.05 


0.17 


Body mass index (kg/m^) 


-0.14 


0.0001 


-0.15 


<0.0001 


Waist-to-hip ratio 


0.02 


0.66 


0.02 


0.63 


Vitamin D intake ^ (lU/day) 


0.13 


0.0003 


0.12 


0.0002 


Calcium intake ^ (mg/day) 


0.15 


<0.0001 


0.13 


0.0005 


Caloric intake (kcal/day) 


-0.10 


0.005 


-0.08 


0.04 


Physical activity (MET-h/week) 


0.20 


<0.0001 


0.20 


<0.0001 


Alcohol (servings/week) 


0.04 


0.34 


0.02 


0.56 


Season at blood sampling, n (%) 


0.40 


<0.0001 


0.40 


<0.0001 


Education, n (%) 


-0.12 


0.001 


-0.13 


0.0004 


Smoking status, n (%) 


0.07 


0.08 


0.07 


0.06 


Total of all colas (serving/week) ^ 






-0.11 


0.004 



' Adjusted for all other determinants included in this table except for total of all colas. ^ Adjusted for all other 
determinants included in this table. ^ Total intakes from food and supplements. " Coded winter = 1, fall and 
spring = 2 and summer = 3. ' Including regular colas and caffeine free colas. 
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Spearman correlations between each determinant of 25(OH)D and plasma 25(OH)D concentration 
are shown in Table 3. The intake of vitamin D and calcium, season at blood sampling and physical 
activity were positively correlated with 25(OH)D concentrations (r = 0.13 {p = 0.0003), r = 0.15 
(p = <0.0001), r = 0.40 (p = <0.0001) and r = 0.20 (p = <0.0001), respectively). Negative correlations 
with 25(OH)D concentrations were observed for BMI, caloric intake and education (r = -0.14 
(p = 0.0001), r = -0.\0(p = 0.005) and r = -0.12 (p = 0.001), respectively). In our data, no correlation 
was observed between age, waist-to-hip ratio, smoking status or intake of alcohol and 25(OH)D 
concentrations. Adding total intake of all colas in the model slightly changed the observed correlation 
between education, calcium or caloric intake with 25(OH)D concentrations (r = -0.13 (p = 0.0004), 
r = 0.13 (p = 0.0005), r = -0.08 (p = 0.04), respectively), but all of the other correlations remained 
essentially unchanged. 

4. Discussion 

In this study, we found that higher intake of colas was associated with lower plasma concentrations 
of 25(OH)D among premenopausal women. Compared to non-consumers, women who drank more 
than three servings per week of colas presented a mean concentration of 25(OH)D that was 12.7% 
lower. In the latter category of colas consumption, 47.8% of women presented 25(OH)D 
concentrations below the recommended 50 nmol/L. Furthermore, 10.9% of these women had 25(OH)D 
concentrations below 30 nmol/L, which are considered as deficient concentrations. In comparison, 
only 16.3% of the women who never drank colas presented concentrations of 25(OH)D below 
50 nmol/L of which only 0.3% {n= 1) were below 30 nmol/L. 

Few studies have evaluated the association between intake of sugar-sweetened beverages and 
concentrations of 25(OH)D [24,25]. Following a two month diet, Garcia-Contreras and colleagues 
observed lower plasma 25(OH)D concentrations among rats consuming colas beverages compared to 
those who drank water [24]. In a cross-sectional study conducted among 411 obese children, Olson and 
colleagues found that both intake of sodas and juices were inversely associated with serum 
concentrations of 25(OH)D [25]. Our observations appear consistent with these studies conducted in 
population other than premenopausal women. One study, among women, has evaluated short-term 
effects of carbonated beverages consumption on serum concentrations of 1,25-dihydroxy vitamin D 
(l,25(OH)2D), the active form of vitamin D metabolized by the kidney and target tissues (such as 
breast, prostate, intestine, skin and bones) from 25(OH)D [38]. Concentrations of l,25(OH)2D were 
lower in the group following a diet high in colas intake (1.4 L/day) compared to women following a 
diet without cola (85.4 nmol/L and 99.8 nmol/L, respectively), but the association was not significant. 
However, the study population was very small (eight women) and the serum concentrations of 
l,25(OH)2D may not be the best biomarker to represent vitamin D concentrations [22]. Some authors 
have brought the idea that the effect of intake of colas on 25(OH)D concentrations was due to a 
replacement of the milk consumption by colas [25]. However, further adjustment for milk intake did 
not materially alter our results, suggesting that this factor does not confound the observed associations 
with 25(OH)D concentrations. 

In our study, we found an inverse association of intake of colas with 25(OH)D concentrations and 
also a non-significant negative correlation between intake of other carbonated beverages and 25(OH)D 
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concentrations. However, the intake of sweet fruit drinks was not associated with 25(OH)D 
concentrations. In comparison to sweet fruit drinks, carbonated beverages contain a higher 
concentration of fructose which is derived from com syrup and used as a sweetener [39]. It has been 
suggested that fructose may have an effect on vitamin D metabolism [26]. In fact, the liver and kidneys, 
the principal sites of vitamin D metabolism, seem adversely affected by chronic fructose intake [40,41]. 
Moreover, one study conducted among rats with compromised renal ftinction has shown that those fed 
with fructose had a reduction of 30%-40% in their 25(OH)D serum concenfrations compared to rats 
fed with glucose [26]. This observation could explain the observed lower 25(OH)D concentrations 
with higher intake of carbonated beverages in our population. Furthermore, we observed no association 
between intake of diet colas and 25(OH)D concentrations (data not shown). However, the statistically 
significant association found with intake of colas was stronger than that of intake of other sweet 
carbonated beverages and thus may not be solely due to fructose. 

A difference between colas and other sugar-sweetened beverages is that they contain caffeine which 
has been shown to increase the risk of osteoporosis [42,43]. To take this into consideration, we have 
made separate analyses of intakes of regular colas and caffeine free colas with 25(OH)D 
concentrations, and found similar negative correlations for both types of colas (r = -0.094, p = 0.012 
for regular colas and r = -0.103,/? = 0.006 for caffeine free colas). Thus, in the present study, caffeine 
does not seem to explain the results. 

In contrast to other sugar-sweetened beverages, colas also contain phosphoric acid [5] and excessive 
exogenous phosphate intake is believed to cause vitamin D and calcium metabolism disorders [27,28]. 
One study has shown that intake of soft drinks containing phosphoric acid induced hyperphosphaturia 
and significant reductions in plasma concentrations of l,25(OH)2D and 25(OH)D in immature rats [44]. 
All of these hypotheses need to be confirmed with additional studies, especially among humans. 

This study has several strengths. The study sample was relatively large and the women were 
recruited in a short period of time. The measurement of plasma 25(OH)D concentrations was done 
with respect to high quality standards and the food frequency questionnaire used to assess intake of 
sugar-sweetened beverages was validated in other studies for its accuracy [33,34]. Several possible 
confounding factors (including vitamin D intake from food and supplements and milk consumption) of 
the association between intake of sugar-sweetened beverages and 25(OH)D concentrations were taken 
into account. Furthermore, the women of this study lived in the same geographic area, thus were 
exposed to the same climatic conditions and the assessments of leisure-time physical activity with the 
season at blood sampling allowed for a good adjustment for sun exposure. Thereby, the difference seen 
in 25(OH)D concentrations was probably not caused by differences in vitamin D intake, milk intake or 
sun exposure. 

This study has also some limitations. Misclassification in the frequency of consumption of 
sugar-sweetened beverages is possible with the use of a food frequency questionnaire. Intake of 
sugar-sweetened beverages over the previous year was an estimation given by each participant. 
However, differential bias is unlikely because answers inaccuracies can hardly be related to 25(OH)D 
concenfrations and if there was misclassification, it could only lead to an underestimation of the effect. 
Also, it was difficult to state on the biological mechanism because the nutrients contained in colas that 
can be responsible for the lower 25(OH)D concenfration (finctose, caffeine and phosphoric acid) were 
already present in several other foods. Despite infroduction of season at blood sampling and 
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leisure-time physical activity in the models, the adjustment for sun exposure was not perfect because 
the information about clothing style, sunscreen use and time spent outdoors other than leisure-time 
physical activity were not available. It is therefore possible that colas drinkers spent less time outdoors 
than non-drinkers and this could lead to residual confounding. Intake of colas may also be viewed as a 
proxy of an unhealthy lifestyle which could include higher BMI, waist-to-hip ratio, intake of alcohol 
and calorie and smoking habit, and lower intake of calcium and vitamin D and levels of physical 
activity. However, in our data, the correlations between such 25(OH)D determinants and 25(OH)D 
concentrations that were statistically significant remained so even after the introduction of total intake 
of colas in the model. Furthermore, this study population was limited to premenopausal women so 
inference to other populations is difficult to make. Finally, the cross-sectional design of this study does 
not allow causal interpretation. 

5. Conclusions 

In conclusion, we found that premenopausal women with higher intake of colas had lower 
circulating concentrations of 25(OH)D. To our knowledge, this is the first report of this association in 
such type of population. Considering the importance of vitamin D concentrations in the maintenance of 
a good health and the possible links of vitamin D deficiency to some diseases, our findings need 
further investigation. In our study population, the intake of colas was not very high, but in other 
populations this consumption can be more significant and may affect 25(OH)D concentrations in a 
greater manner. Considering the high consumption of this type of beverage in the general population, 
further studies are needed to elucidate the findings of the present study. 
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